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Abstract. The solutions to a large class of semi-linear parabolic PDEs 
are given in terms of expectations of suitable functionals of a tree of 
branching particles. A sufficient, and in some cases necessary, condi- 
tion is given for the integrability of the stochastic representation, using a 
companion scalar PDE. 

In cases where the representation fails to be integrable a sequence of 
pruned trees is constructed, producing a approximate stochastic repre- 
sentations that in some cases converge, globally in time, to the solution 
of the original PDE. 



1. Introduction 

This paper considers stochastic representations for solutions to a large 
class of semi-linear parabolic PDEs, or systems of PDEs, of the type 

(1.1) d t u = Au + ^(u)+f, 

where A is an operator with a complete set of eigenfunctions, & is a polyno- 
mial nonlinearity in u and its derivatives, and / is a given driving function. 

In short, the solution u is expanded into a Fourier series using the eigen- 
functions of A. This yields (as in spectral Galerkin methods) a system of 
countably many coupled ODEs for the Fourier coefficients. This ODE sys- 
tem is then solved in a weighted £°°-space, via an expectation over a tree 
of branching particles. The rules for the branching and dying probabili- 
ties arise from the particular PDE being studied. Moreover the PDE deter- 
mines an evaluation operator R t , acting on the tree of particles rooted at 
each Fourier mode k, so that, under integrability assumptions, the (suitably 
weighted) k'th Fourier mode %k(0 is given by 

Xk(0 = Eft(<r k )]. 

This is precisely the method of Le Jan and Sznitman [ 8 1 (later extended in 
Bhattacharya et al. |3 1, Waymire ifTTll . or Ossiander lfT3lO where they treated 
the Navier-Stokes equations in R 3 . Earlier papers connecting branching 
particle systems to PDEs (for instance Skorokhod 1T31 orlkeda, Nagasawa, 
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and Watanabe |7 1, and later McKean lITTlO use branching coupled with a dif- 
fusion, and the stochastic representation is derived directly without Fourier 
series, so that the linear operator A is limited to generators of diffusions. 

Our first aim is to show that this method applies to a large range of equa- 
tions. In Section El we present three representative examples with quadratic 
non-linearities. We comment on further generalisations in Section 1231 Ba- 
sically any system of parabolic PDEs with polynomial nonlinearities in the 
derivatives is admissible. 

One major drawback of the stochastic representation is that it often fails 
to exist for large times t, although the solution to the PDE may still exist. 
The problem is that Rt(Tk) may fail to be an integrable random variable 
for t > to. When deriving a system of ODEs in £°° space, there is consid- 
erable freedom in the choice of weights for the Fourier coefficient under 
which integrability can be established. See Bhattacharya et al. [3j for an 
extensive discussion in the case of 3Z)-Navier Stokes. One way to check in- 
tegrability is to establish a scalar comparison equation. The finiteness of the 
comparison equation implies the integrability needed for the stochastic rep- 
resentation to hold. The comparison equation is independent of the weights 
and represents a worst case scenario with super-linear (explosive) growth. 
It typically completely ignores most of the structure of the non-linearity in 
the original PDE. 

In Section|3]we establish the stochastic representation under the assump- 
tion that it is integrable. In Section |4] we investigate the comparison equa- 
tion. This typically shows the representation is integrable at small times, 
or, when there is no linear instability, for all times with small data. For 
some classes of equations, for example U-Burgers equation, we obtain a 
necessary and sufficient condition for the integrability of the stochastic rep- 
resentation, independent of the choice of weights in Fourier space. 

Our second aim is to present an approach to treat cases where integrabil- 
ity fails. The key point is to get rid of the smallness condition on the initial 
data and the forcing, in order to find a stochastic representation that is global 
in time. In Bhattacharya et al. [3 1 the branching trees are pruned after n gen- 
erations. This gives a stochastic representation of a Picard iteration scheme 
converging to the original PDE, but, as stated in [3|, the existence of the 
expectation is equivalent to the convergence of the Picard iteration scheme. 
In another approach, Morandin lfl2l suggested a clever re- summation of the 
expectation in order to improve the convergence for large times, but he was 
only able to rigorously verify the global convergence of his method in a 
simple example where (11.11) is a one-dimensional ODE. Our approach is to 
construct sets Q. n , with P[Q.„] | 1 so that 

Xk(0= limE[^(T k )la n ]- 

This treats the expectation somewhat as a singular integral, where we have 
to be careful how to cut out the singularity. 
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The method we use, explained in Section is to construct a pruned 
branching tree which will agree with Tt on Q. n . The expectation for 
the pruned tree 

E[^(^ k ( " ) )] = E[^(^)lnJ 

is well defined and represents the k th Fourier mode of the solution to a 
semi-implicit approximation scheme of the type 

d t u (n) =Au (n) + ^(u (n \u in - l) )+f. 

We then use PDE techniques to verify that the approximation scheme con- 
verges to a solution of the original PDE. Although there are general results 
for the convergence of such approximations (cf. for example Bj0rhus and 
Stuart [4J) the assumptions are usually quite restrictive. Since stronger ar- 
guments are model specific, we present the arguments only in two special 
cases, namely for a simple quadratic ODE and for Burgers equation. We 
believe these examples illustrate that the method potentially dramatically 
extends the range of PDEs for which there is a global stochastic represen- 
tation. A global result is essential if one wants to study a stochastic rep- 
resentation of the long time behaviour of solutions, for example in terms 
of stationary solutions or pull-back fixed points. Only in a very simple 
framework of small initial conditions and uniformly small forcing is it cur- 
rently possible to derive such results (for the 3D Navier-Stokes equations 
see Bakhtin [2J and Waymire lfT71 ). The extension of such results to non- 
trivial cases and the relation with the pruned representation are the subject 
of work in progress. 



2. Abstract setting and examples 

We first present an infinite system of ODEs involving a quadratic non- 
linearity. The system is indexed over k G Zf We then discuss several 
examples of PDEs on the torus [0,27t) rf and recast their Fourier transforms 
into our abstract ODE setting. We do not present the highest generality 
possible, but focus instead only an equation with one quadratic nonlinearity, 
one linear instability and one forcing term. We comment in Subsection 
12.51 on a large number of possible extensions, including other domains and 
boundary conditions, multiple forcing terms and additional nonlinearities, 
possibly of higher order. 

2.1. The general system of ODEs. We consider solutions %{t) :Z d ^C 
to the following infinite dimensional system of C'-valued ODEs 

(2.1) % k = ?l k -Xk + Cfp k X k + C b £ %,l,m#k,l,m(Xl,Xm)+dkYk 

l,meZ d 

with k E Z d . The constants > (which will determine the rate of particle 
evolution), p^, gk,i,m> £ [0, 1] (which will determine the probabilities of 
flipping, branching and dying), and Cf, C\, > (the flipping and branching 
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constants) are fixed, as are bilinear operators 5k,i,m : C xC->C satisfy- 
ing 

|£k,l,m(X,X')l < Ixllxl 
for all x,x' £ C r . The choice of these constants will arise from the Fourier 
transform of the PDE being studied. We assume throughout that 

(2.2) Pk + qk + d k =l forallkeZ rf , 
and 

(2.3) pk ^ 0, qk -> 0, as |k| -> °o 
where 

The data for the equations consists of a time dependent forcing y = {7k (0 : 
k G Z d , t>0} and an initial condition x(0) = {Xk(0) : k G Z d }. We consider 
the above system in its mild formulation, that is for given data we look for 
measurable 1 1— > Xk(0 £ C satisfying, for k G Z d , 

(2.4) Xk(0 = e-^Xk(O) + J\ k e- K{t - s) [c f pMk(s)+ 

+ C b <lk,l,mBk,l,m(Xl(s),Xm(s))+dkyk(s) ds. 

Note that we need some regularity of Xk> in order to make (12.41) well defined. 

Remark 2.1. There is considerable flexibility when choosing the constants 
in the ODE system (12.11) . For example, we can adjust the probabilities pk, 
<7k,i,m 5 and dk by adjusting the constants Q,, C/ and considering modified 
forcing data 7. In particular, in an equation where the probabilities do not 
add up to 1 in (12.21) . it is always possible to adjust dk and the forcing data 
so that this constraint holds. Similarly, an equation with Cf and Q, replaced 
by bounded functions of k can be recast into the form (12.1ft by forcing the 
k dependence into the probabilities pk, qk, an d dk- 

2.2. The d-dimensional Burgers equations. Consider solutions u(t,x) G 
R rf , for t > and x G [0, 2n) d to the Burgers system 



(2.5) 



d t u — Au + (u ■ V)u = f, 

U(0) = M°, 

with periodic boundary conditions, where / is an external forcing. We re- 
strict ourselves to periodic boundary conditions, as the nonlinearity is easy 
to compute in the Fourier basis. Nevertheless, other kinds of boundary con- 
ditions, for instance like Dirichlet or Neumann, can be treated in a similar 
fashion (cf. section l2~5l) . 
If we expand the solution 

u(t,x) = £ Uk(t)e ik -\ 
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the equation reads in the Fourier coefficients as 

u k = -|k| 2 w k -i (>i-m)wm + A, 

l+m=k 

The sum is over all 1, m G Z d satisfying 1 + m = k. Define a weight function 
Wk = 1 V |k| a , where a > will be chosen shortly, and set Xk = Wk"k- Then 



(2.6) 




|k| 2 Xk - ili +m=k S(» • ^)% m +/ k w k , 
"k(0)w k 



Note that the mode uq has no linear dissipation. Below we will add and 
subtract XqXq to the equation for %q, which introduces a linear instability 
but which allows us to write the equation in our desired abstract form. We 
note that in dimension d = 1 this trick is unnecessary: the equations for the 
zero* mode decouples, in that it simplifies to uq — fo, and it is then possible 
to reduce the problem to the case fo = uq = 0. 

We now show one way to recast (12.61) into the abstract form (I2.ll) . For 
given Cf, Cj,, > we define 




k^O, 
CJ 1 k = 0, 



whenever 1 + m = k (and zero otherwise). Lemma I2T21 below ensures, pro- 
vided we choose a > max{^-, d — 1}, that q k = Ei m <7k,l,m < °° an d that 
q k — ► as |k| — > oo. Thus by taking Q,, Cy sufficiently large we have 
that /?k + <?k < 1 an d it remains only to define d k = 1 — /?k — <7k and Yk = 
(AwkAk^k) f° r k G Z rf . Again, there is considerable flexibility in these 
choices. 

Lemma 2.2. For all a, y > wiY/i a + y > d there exists C = C(a,y) < °° 
so that, for all k G Z , with k 7^ 0, 



'c(l + |k|) p , ifa^dandy^d, 
m\ a \\\y ~ 1c(l + |k|)~ p log(l + |k|), ifa = dory=d, 



< 



l+m=k 

where [3 = min{oc, y, a + y — d} and the sum is over all indices 1, m in Z d 
satisfying the given constraints. 

One way to prove this lemma, whose proof is omitted, is to compare 
above and below by suitable continuous integrals. 

2.3. Two dimensional Navier-Stokes equations. We briefly treat the two 
dimensional Navier-Stokes in its vorticity formulation, since this will be 
used in section|4|as an example where the comparison equation yields exact 
statements about the integrability of the stochastic representation. 
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In dimension d = 2 the vorticity £ = curlw is a scalar and satisfies, on the 
torus [0,2ft) 2 and with periodic boundary conditions, 



(2.7) 



3£-as+(«.v)$ = /, 

«0) = 



where u is the solution to the Navier-Stokes equations. The Fourier coeffi- 
cients satisfy the following system, 

4k = -|k| 2 ^ k + £ 

l+m=k I 1 ! 

where 1^ = (fe, —h). For simplicity we shall assume that /o = and the 
vorticity has mean value zer o and is omitted from the system. 
We then set %k = |k| a ^k for some a > j. For Q, > we then define 



k • 1 

l k = |k| 2 , 5k,l,m (%,%') = j^pjXX', 4k,I,m = C b l 



|k| 


a-2 


k 


i 1 ! 




1 


a+2 


m 


a 



for all k, 1, m 6 Z 2 satisfying k ■ l- 1 ^ and 1 + m = k (and zero otherwise). 
Lemma l2~2l ensures that <?k < 00 an d that — > as |k| — ► °°. Taking Q, 
large enough we have that % < 1 (note that here we may take = 0). So 
the recasting is complete if we define = (|k| a_2 /Jk)/k- 

2.4. A surface growth equation. The final example illustrates the change 
in weights needed for a higher order equation and the need to consider linear 
instabilities. In particular, the linear operator does not generate a diffusion. 
Therefore, the Fourier transform is necessary for the stochastic represen- 
tation. Consider the following scalar equation arising in some models for 
surface growth, 

dtu — — aiA 2 w — a2Au — asAlVul 2 + a,4\Vu\ 2 + f , 

with periodic boundary conditions on [0,2n) d , with d = 1,2 and a ; > 
for i = 1,2,3. See Raible et al. lT4l for the derivation of the model, and 
Blomker et al. [5| for a rigorous mathematical treatment using PDE tech- 
niques. For simplicity, we assume a$ = and that the mean value / u(t,x)dx 
is zero, allowing us to omit the coefficient uq. 

The equation for the Fourier coefficients is given by 

(2.8) iik = — «i |k| 4 Wk + «2 |k| 2 Wk + <33 |k| 2 £ (1 ' m ) M i M m + A- 

l+m=k 

We set %k = |k| a Wk for a > 0, with a > maxjd, 1 +d/2}, and then choose, 
for all k 7^0, 



and 



K = a,\k\\ p k = -C7 1 |kr 2 

a\ J 



.. 1-m , _ 1 «3 Ikl 01 2 |1 - ml 

5k,l,m(X,X) = Tr^nXX 5 4k,I,m = <V 



].mr" J " ai|l| a |m| a 
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with 5^ i m and <?k,i,m equal to zero if 1 • m = or 1 + m ^ k. Lemma I2T21 
guarantees that Pk + Qk< 1 when Q,,C/ are taken large enough and we can 

|k|«- 4 

define = 1 — — /?k and Yk = ' ' , jk to obtain a system in the form 

" 1 "k 

dTTT) . 

2.5. Discussion of extensions and generalisations. We now list a number 
of possible extensions to the our basic system (12.11) for which modified tree 
representations will hold. 

2.5.1. PDEs in general domains with other boundary conditions. If there 
is a complete countable set of L 2 -eigenfunctions (e k ) ke j^ °f A in which to 
expand solutions as u — YX=\ u k e k, one can recast PDEs in general domains 
with various boundary conditions into a suitable ODE setting. This is simi- 
lar to spectral Galerkin methods. Consider for instance 

d t u = Au + B(u,u) 

for a bilinear operator B, and suppose that Ae k = X^k- Then 

d t u k = X k u k + ( B (e m ,ei),e k ) L 2U m ui. 

m,l=l 

This can easily be transformed into the general system (12.11) by choosing 
appropriate weights. This would cover our earlier examples, Burgers equa- 
tion, Navier-Stokes or the surface growth equation, in a regular domain 
with, for instance, Dirichlet or Neumann boundary conditions. Note that 
(12.11) is now an R- valued system posed in £ TC (R) which is indexed over N 
instead of Z. 

2.5.2. Polynomial non-linearities. More general polynomial non lineari- 
ties, or several non-linearities, lead to branching systems where particles 
split into a larger number of descendants. Even analytic non-linearities can 
be handled, with the absolute values of the power series coefficients con- 
trolling the branching probabilities. Note also that a general first order term 
of the form Li,Pk,i#k,i(%k), for linear 5k ; i : C — > C, can be thought of as a 
branching event with a single offspring. This kind of term arises, for exam- 
ple, when the original PDE contains a multiplication operator fu for a 
fixed function /. 

2.5.3. Multiplicative forcing. A non-linear forcing term F(u, f), again with 
polynomial F, can also be recast into a branching system of ODEs. This 
leads, say in the quadratic case, to time dependent bilinear operators 5k,i,m 
whose values depend on the forcing 7k(0> i- e - we obtain terms like a sum 
over <7k,i,m5 k ,i,m(Yi,Xm) in equation (O)- 
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3. The branching particle representation formula 

3.1. Existence and uniqueness. The next theorem shows that there is a 
unique local solution to (12.11) taking values in the space t°(C) of bounded 
families («k)kez rf °f elements of C, with the norm ||a||oo = sup keZ </ \a\i\, 
with | at | = y/ak • a^ the norm in C. We give a simple deterministic proof, 
but there is also a more probabilistic proof available (see Corollary 14.31) . in 
the spirit of Le Jan and Sznitman [ 8 ] . 

Theorem 3.1 (Unique local existence). Assume that 

Z(0)6f(C r ), YeL oo ([0,r],r(C r )) farallT>0. 

Then there exists a time Tq > 0, depending only on %(0), y, and the constants 
appearing in the equation, such that the mild formulation (12.41) has a unique 

solutionxeL? oc ([0,To),(C r ) Zd )- 

Moreover, we have either Tq = °° or \\%(t) ||«, — > °° as t — > Tq. Finally, if 
the functions t i— > Jk(t) are C k , then t i— > %k(0 are C k+l in time and solve 
equation (12. II) . 

Proof. The proof is a rather standard application of the Banach fixed point 
theorem. Let B be a ball of radius R > centred at the constant function 
with value %(0), in the space L oo ([0,^],£ oo (C'")). For % e B define F(%) by 
the right-hand side of Q>. Then for Rq=R + ||%(0) ||oc, 

l^(x) k (0-Xk(o)| < (llxCO)!!^^^^^^^^^ ||y(0l|oo)(l -e-^*). 

If we choose R > ||%(0) ||„ + sup \\y(t) ||oo and ^ small enough we see that F 
maps B into itself. Here we have used assumption (12.31) to control the large 
|k|s. Moreover, if y} and % 2 are in B, then for t < U, 

|[^(% 1 )-^(% 2 )]k(0l<(QPk + 2C^ %)(l-e" V *)sup||% 1 (0-% 2 (0l|cc. 

Hence i 7 is a strict contraction in L co ([0,t i ],£°°(C r )) if we choose ^ small 
enough. Here we need again, for large |k|, the assumption (12.31) . 

The assertion for the time Tq follows in a standard manner by gluing 
together local solutions. The continuity of t — > %k(t) follows from the 
mild form (I2.4I) . It is even differentiable with bounded derivative. The C - 
regularity follows by differentiating (12.41) and the higher regularity follows 
from differentiating (12. 1 1) . □ 

A simple global existence result can be proved under the assumptions of 
linear stability and small data. 

Proposition 3.2 (Global existence for small data). Under the assumptions 
ofTheorem \3.1\ assume that there exists 8 > so that 

dklYkl < 8(1 -C f p k ) -C h d 2 q kl for alike Z d . 

Then for each initial condition \\%(0) ||oo < 8, there is a global solution % to 
equations (12.11) satisfying sup f>0 ||%(?)||oo < 8. 
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Proof. Let v k (r) = |%k(0 1 • When v k (t) ^ and v x {t) < 8 for all X G Z d one 
has the estimate 

(3.1) 3 f v k < -^k(l-C/p k )v k + ?i k (C fc <5r k 8 2 + Jk|yk|)- 

The assumption implies that the right hand side of (13.11) is negative when 
v k = 8 and global existence follows from a comparison argument for one- 
dimensional ODEs. □ 

3.2. The branching tree. We now give a construction of the branching 
process that will be used to represent the solutions of d2.ll) . We will label 
particles of the process with labels taken from the set J* = \J^ =0 {0, 1,2}". 
The history of a particle a = (oci, . . . , a M ) can be read off by interpreting 
aj = as the flip at generation j, and a 7 = 1 (or 2) as being child 1 (or 2) 
in a binary branching event at generation j. 

For a G {0, 1,2}" we write |a| = n which we call the length of the label. 
We write a = for the single label of length zero. When a = (ai , . . . , a„) 
we write a|j for the label ct\j = (oci, .. . ,0C/) of its ancestor at genera- 
tion j G {0,1, 1} (and set a|0 = 0). For i G {0,1,2} we write 
(/, a) for the label (?', cci , . . . , a n ) and (a, i) for the label (oci , . . . , a n , i) (or 
(/, a) = (a,/) = (/) if a = 0). We construct the branching particle sys- 



FlGURE 1 . A tree with branches, deaths (•) and flips (o). 

terns on a probability space equipped with the following independent fam- 
ilies of I.I.D. variables: (E a ) ae: y exponential mean one variables (that 
will control the overall rates of branching and flipping); (£/ a ) ae ,y uni- 
form [0, 1] variables (that will control whether a particle flips, branches or 

dies) and ((^(k),!^)) 

aej^.kez d random variables with distribution 

i 5 [y'a 1 ' ) (k) = 1, Yap(k) = m] = g ki i im (which will control the positions of 
the two offspring of a particle that branches). 

We now define a system (K a ,x^,x^) a(Ei y of particle positions, birth and 
death times, inductively over the length n = |oc| of the labels. Fix k G Z rf 
and set Kq = k, = and x^ = XT E^. Assume that the positions, birth 
and death times have been defined for |oc| < n. Then, for a of length n + 1, 
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define birth and death times 



B _ D D_ B , p 

x a _ x a|rc' x a — x a + A K a |„ a ' 



and the particle positions 

[ K a i„, a n+ i = 0, 

Ka = \ ia (K a | n ), CC n+ i = 1, 

l^a ^(K a | n ), a n+ i = 2. 

This defines a complete tree of all possible branching and flipping particles 
rooted at k. In the desired evolution the particles will choose whether to 
flip, branch or die according to the probabilities p^, q^, d^. 



Figure 2. The construction of the tree. At each event time 
x, there is a random selection between either death, flip or 
branch of two new particles to states 1 and m (depending on 
the state k of the parent particle) 

We now define indicator variables (Ia)ae,y to decide whether a particular 
branch has survived. Define 4 = 1 and, for a of length n + 1, 

if a„+i = and U a E [0,/^ ], 
/«= <( 1 if oc+i G {1,2} andU a e [1-^,1], 
otherwise. 

Fix an isolated cemetery state A and define, for jocj = n, 



Ka ifrr)=i/au = i, 

A otherwise. 



The collection Tk = (K a , x„, i^) ae j? now defines our branching tree rooted 
at k. It lives in the space defined by 

^=((Z rf U{A})x[0,oo)x[0,oo)y / . 

We denote the law of Tk on by Pk. 

The descendants of any one particle in the tree form a new tree. To make 
this precise we define shift maps 7t, • : & — > & , for i = 0, 1,2 as follows. For 



A BRANCHING PROCESS REPRESENTATION. 



11 



1 = (k a ,5 a ,? a ) aei y G £F we define a new tree 7t,-(T) by 

7l i'(' r ) = (k(i,a)» s (i,a) -h,t{i,a) -h)aeS- 

The tree TC ; (T ) is meant to be the tree of descendants of the particle labelled 
(i), with their birth and death times shifted so that particle (/) is born at 
time t = 0. The construction of the branching particle system from I.I.D. 
families implies the following lemma. 

Lemma 3.3. Let = (K a ,x^,x^) ae j? have law Pk. Then 

1. conditional on {x^ G ds, K/q) = k} the tree 7To('7k) ^fls dze law Pk/ 

2. conditional on {x^ G ds, Km = m,K( 2 ) = 1} ^ e frees Tti(Tk) ^d 
7l2(Tk) are independent and have laws P m and Pi. 

We want to ensure that the tree has only finitely many branches before 
time t. Define iV[ 0)f ] : & -> N by iV[ 0i? ](T) = |{a G : s a < 01' mat is me 
cardinality of the set of particles born before time t. 

Lemma 3.4. Under Pk the variables Mo,*] are almost surely finite for all 
f >0. 



Proof. Let Pb(t) = Pk[N[Qj] < °°]- By conditioning on the values of Xj 
Km), K( 2 ) and using Lemma 1331 




Hence, (Pk(t) : k G Z d , f > 0) is a bounded, real- valued solution to the equa- 
tion (12.11) with forcing y = 1 and bilinear operators Bk\ m (%,%') = By 
Theorem 13. 11 there is only one solution, namely Pk(t) = 1 for all k, t. □ 



A simple criterion that ensures that the branching process becomes ex- 
tinct with probability one, that is K a = A for all large |oc|, is that 

(3.2) a k < d k and p k < 1 for all k G Z d . 

Indeed the number of particles alive at time t is an integer valued pro- 
cess whose successive values, under the condition (13.21) . form a sub-critical 
branching process. Therefore it eventually reaches zero. The number of 
values k G Z d taken by particles before this extinction is almost surely fi- 
nite. The conditions that X k > and pk < 1 ensure that the extinction time 
for the branching particle system is almost surely finite. Note that, as ex- 
plained in remark |2~T1 we can always choose the system (12.11) in such a way 
that (Q holds. 

3.3. The evaluation along the tree. We now fix a forcing functions y and 
an initial condition %(0). We wish to define evaluation maps R t : & — > C 
for t > 0, which will depend on y and %(0). These will satisfy a recursive 
property that allows them to be calculated backwards along the tree. 
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For the sake of simplicity, we introduce the following abbreviations: 
given a branching tree T = (k a ,%,r a ) ae j?r and a particle labelled ae/, 
with k a 7^ A, we say that the particle has a 

death: if k (Oj0) = k ((X)1) = k (oC;2) = A, 

flip: if k (a ) # A and k (0)1 ), k (oC;2) = A, 

branch: if k (ot;0) = A and k (ct;1) , k (0j2 ) 7^ A. 
Under each probability P k , every particle a for which K a 7^ A must do 
exactly one of the above three possibilities. 

Lemma 3.5. There exist a family of maps R t : — > C, /or ? > 0, satisfying, 

when N\q r ] (T ) < <*>, implicit formula 

(3.3) 

Xke(0) ?0 > 

Yk (f — t<D <t, death at 0, 

C / i? f _ f0 (7U O ('r)) t & <t,flipat<b, 
k C fc fi k0;k(1);k(2) (i? f _ f0 (7Ui(T )),^_, (7U 2 (T ))) < f, oranc/i a? 0. 

/Voo/ Informally, since the tree is finite when M 0;f i < 00 the value of R t (T ) 
can be calculated backwards along the tree, starting at time s = t and work- 
ing back to time 5 = 0: evaluate the initial condition %(0) at any particles 
that are alive at time t, evaluate the forcing function y(s) at any particle that 
dies at time s < t, and apply the bilinear operators at the times of branching 
events. 

For a careful proof one can define a sequence of approximations R^ p t prx) in 
the following way: 

^apprx, (<r ) = / Yk (?- fo) if < t, death at 0, 
I 1 otherwise, 

and^ p ; x) ? (T) is given by 
(3.4) 

'Wo) t<b>t, 

Yk (t -t<b) t®<t, death at 

CfB%R(M*)) ?0<?,fli P at0 
C,5 k0 , k(1) , k(2) (J^(JCi(^ )),«Sq,(3C2(^r ))) < f, branch at0. 

If Af[o. f ] < 00 then only finitely many iterations are needed and Rt(T) = 
Rn7 x) (T) for all large n. □ 

In some cases the evaluation can be written more explicitly. Let F(t) 
(respectively B(t)) be the number of particles that have flipped (respectively 
branched) before time t. Let D(t) be the set of labels of particles that have 
died strictly before time t. 

Consider the special case where r = 1 and that all the bilinear forms B^ \ m 
coincide with the usual product in C. Then the evaluation is given, almost 
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surely under P^, by 

(3.5) R t (T)=C B b {t) C F / t) J] W-'a) EI Xk a (0). 

aeD(t) a:te[s a ,t a ) 

In the general case, we can only verify, under P^, that 

(3.6) \R t (T)\ <<*®<f® J] II 1^(0)1, 

aeD(t) a:te[s a ,t a ) 

and that equality holds in (13.61) if |5k,i,m(%,%')l = Ixllx'l f° r all k,l,m and 
Xi X • 

3.4. The representation formula. Consider an initial condition %(0) e 
f°(C), and a forcing y e L°°([0, r],f°(C)). The representation formula 
for solutions of d2.ll) . when the expectation exists, is given by 

(3.7) Xk(0 = Ekfc], kGZ rf . 

Theorem 3.6. Suppose that there exists C = C(y,%(0), T) < °° so that 
Ekl^fl <C for alike Z d and all t e [0,T]. 

Then % defined in (13.71) w unique L°°([0, r],£°°(C'')) solution of problem 
(12.11) for data y, %(0). 

Proof. Note that uniqueness follows from Theorem l3.ll Fix a k G Z rf . Con- 
ditioning on the values of To, Km), Km, Km) and using lemma 1331 leads 
immediately to the mild form of the equation (I2.4I) . The uniform (over k) 
integrability is necessary to show that the sum over 1, m converges. □ 

In the next two sections we discuss how to check the integrability as- 
sumption and what to do if it fails. We also see, what happens if the solution 
fails to be in 



4. The comparison equation 

4.1. The comparison equation. The comparison equation for system (12. II) 
is formed by taking the norm of the data |Xk(0) | and |yk| as new data for the 
system 



(4.1) 



Xk = ^k [-Xk + C/^kXk + Cfc^ ^z^k^mMm + ^klYkl] 

Xk(0) = |Xk(0)| forkeZ^. 



We now look for non-negative real solutions Xk(t). 

We also define a modified evaluation operator R t on 2? by the implicit 
formula (13.31) where we use the new data |%k(0)| and |yk| and the bilinear 
operators are replaced by #k,l,m(X)X0 = XX'> the normal product of real 
numbers. Then R~t(T) > and formally we expect that 

(4.2) % k (t) = E k [R t ], keZ d 

should solve the comparison equation. 
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The next theorem confirms this and shows that a finite solution to the 
comparison equation (14.11) is a sufficient, and sometimes necessary, condi- 
tion for the tree expectations Ek[i? ? ] to exist. 

Theorem 4.1. If the expectations in (14.21) are finite for all t E [0, T] and 
k E Z d , then they define a mild solution to the comparison equation (14.11) 
for which t — > Xk(0 is continuous on [0, T}. 

Conversely if there exists a finite mild solution of (14. II) . that is Xk(0 < 00 
for t E [0, T] and k E Z d , then the expectations in (14.21) are finite for t E 
[0, T\, and they define the smallest positive solution of (14. II) . 

Finally, the comparison E^[\R t \] < Ek[/? f ] holds, with equality whenever 
|#k,i,m(X, X') I = Ixl |X1 for all k, 1, m and %, 

Proof. For the first claim of the theorem, condition on the values of Tg, K( ) , 
Km, K( 2 ) and apply Lemma l3~3l to see that the expectations %k(t) = Ek [Rt] 
satisfy the mild form of the comparison equation. Moreover the mild form 
of the equation shows that e^ kt Xk(0 is continuous and increasing in t. Note 
that the convergence of the series in the mild formulation is not a problem 
here, because due to positivity, we can use monotone convergence. 

For the second part of the theorem, let % be a mild solution of the com- 
parison equation (14.11) in [0, T] with data |Xk(0) | and |yk | . Define a sequence 
of evaluations on the trees as follows: set /?Q X f pl) (T , %) = Xk (O and for each 
n > 0, 
(4.3) 

'I W0) | ? 0>^ 
|7kfl (t-tn)) | <t, death at 0, 



C / 47\(7I ('T),X) k><t, flip at 0, 

la^7\(7 Il (^),x)^7\( 7I 2('T),X) f<t, branch at 0. 



(In the language of next section, the evaluation correspond to a prun- 
ing of the tree after n generations and the expectation E^R l ^f' ) (T ,%) will 
solve a Picard iteration scheme for (14.11) ). 

Note that, upon dying, flipping or branching, particles of length n are 
evaluated using the true solution %. Inductively one checks, by conditioning 
on the first event, that for all n > 

(4.4) E k [^7"(<r,x)] =&(t). 

Since N[ ^ < °° under P k we have that (T ) R t (T) almost surely. 
By Fatou's lemma and (l4~4l) we find that E k [^ r ] < Xk(0 < °°- 

The third claim of the theorem is immediate from the upper bound (13.61) 
and the fact that it is an equality under the conditions given. □ 

Remark 4.2. Note that in the above theorem, and its corollary below, we do 
not insist the solutions are bounded in t° . 

The first two parts of the above theorem show that, when there exists a 
finite mild solution % to (14.11) . the function defined by Ek[/? ? ] is the smallest 
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solution to (14.11) lying below %. Note in the case of t° solutions there is 
uniqueness of solutions, as in Theorem l3.ll 

As in Le Jan and Sznitman JSJl, it is possible, when there exists a finite 
mild solution % to (14. II) . to show that n — > R i °J ,) (f ,%) is a non-negative mar- 
tingale (with respect to a natural filtration along generations of the tree). 
Uniform integrability of this martingale would then imply that %k(t) = 
E k [R t }. 

Corollary 4.3. Under the conditions of either the first or the second part of 
theorem Wl\ the expectations %k(0 = E k [/?r] are well defined for t G [0,T] 
and k G Z d and form a mild solution to (I2.ll) . Moreover, such a solution is 
unique among all mild solutions %' such that 

(4.5) |j£(f)| < E k [^], forallkeZ^te&T}. 

Proof. The expectations E k [/?j] are well defined by theorem |4~T1 as \R t \ < 
R t . By conditioning on the first event as before they will solve the mild 
equation. Note that in this case the convergence of the sums in the mild 
equation is guaranteed by the finiteness of the comparison equation. 

Let %' be a mild solution verifying (14.51) and define a sequence of eval- 
uations R { ^\T ,%') for n G N as in the proof of previous theorem, that is 
R { °J\<T ) = % kf) (r) and, for all n > 1, Rff{T ) is defined as in formula (IQl 
with data %' and y and with products 5k,k (1 ),k( 2 ) m the place of usual product. 
By assumption (14.51) and an argument similar to (13.61) it follows that 

\Rfp{T^)\<R^\<T,t), 

where % k (0 = E k [# f ] and ^f Pl) (T,%) are taken from the proof of Theo- 
rem |OJ Moreover, as in that proof, we can show inductively that % k (?) = 

E k [/?g pl) (T,x')]- 

We next note that R%f ] (r,x') =R t (T) and R^Ti? »5c) =^(' r ) on the 
set ft v = {#[<>,,] (<T) < n}. Thus, 

Ek[^(<r,X)losj = % k (0 - E k [^7"(T ,x)la„, t ] 

= E k [^]-Ek[^%J = Ek[^l^ v ], 

and therefore 

|xL(0 - E k [^]| < E k K7'>(T,xO -*| = Ektl^T'^^O -^IHJ, 

< E k [(47"(<r ,%) +R t )l&J = 2E k [R t l a <J. 

Letting n^»we conclude that %' = %, the solution given by the proba- 
bilistic representation. □ 

4.2. Examples. We can remove the weights used to cast the equation into 
our abstract form and rewrite the comparison equation as equations for the 
Fourier coefficients of a scalar PDE. 
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Consider the Burgers equation example discussed in section 1231 Defin- 
ing i<k = w k *%k we obtain a comparison equation of the form 

u k = -\k\ 2 u k + H"l"m+|/k|- 
l+m=k 

"k(0) = \u^\ 

which in the space coordinates corresponds to the scalar equation 

d t u = Au + u( — A)?u + f 

where / has Fourier coefficients \fk\. Note that this scalar comparison equa- 
tion is independent of the choice of weights (called majorizing kernels in 
Bhattacharya et al. Q). 

For the two-dimensional Navier Stokes equation discussed in section 1231 
the comparison equation for ^k = |k|~ a %k takes the form 

Ck=-|k| 2 lk+ £ l ^U m +\A\, 

l+m=k I 1 ! 

which does not have a nice expression in the space variables. 

For the surface equations discussed in section IZ4l the comparison equa- 
tion becomes 

d t u = —cnA 2 u — a2Au — A\(—A)?u\ 2 +f, 

where the forcing / has Fourier coefficients |/k| . 

Whenever there is a solution to these scalar comparison equations with 
finite Fourier coefficients we obtain the existence of mild solutions to the 
corresponding abstract ODEs given by the stochastic representation (13.71) 
This in turn is equivalent to the existence of solutions to the original PDEs 
with finite Fourier coefficients. 

However all three scalar comparison equations have quadratic growth and 
it is possible to show, for example for zero forcing and large enough initial 
data, that the solutions explode in finite time. See for example ffOl and the 
references therein for the case of branching with diffusion. 

Remark 4.4. In the case of the 2d Navier Stokes, the Id Burgers, or the 
surface equation, the equality in the last part of Theorem 14 . 1 1 holds . This 
implies that the stochastic representation Ek[/? f ] is well defined as the ex- 
pectation of an integrable variable, if and only if the corresponding compar- 
ison equation has a solution with finite Fourier coefficients. In particular, 
for any suitable weight, the representation will fail to exist at the same time, 
once a Fourier mode in the comparison equation becomes infinite for all 
solutions. 

5. The pruned approximation 

5.1. An ODE example of the approximation scheme. We first explain the 
main ideas of the approximation scheme on a simple example, namely the 
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equation u = — u + u 2 . The solution can be given by the stochastic repre- 
sentation u[t) = E [^(O)^'], where N t is the number of particles at time t of a 
simple rate one branching process starting from a single particle at time 0. 
It's easy to verify that the representation is well defined for all time t > if 
and only if |w(0)| < 1, in that the variable 1^(0)1^' becomes non-integrable 
for large t when |w(0) | > 1, while the solutions of the equation blow up only 
if«(0)>l. 

We now give a modification of the branching process. Give each particle 
a label from the integers N. Particles still branch at rate 1 but a particle with 
label n produces two offspring, one with label n and one with label n — 1 . 
When a particle of type tries to branch it simply dies. Start with a single 
particle with label n and let N t (n) denote the number of particles at time 
t. Set u n (t) = E[w(0) w 'M] for n > and w_i = 0. Then u n {t) solves the 
following semi-implicit iterative scheme 

u n = —u n + u n -iu n , a„(0)=w(0), forn>0. 

It is straightforward to check that u n (t) is well defined for all n and t. More- 
over, u n converges to the solution u(t) of the original problem for each initial 
condition u(0) < 1. This yields the stochastic representation 

u(t) = lim E \u(0) N '^] 
valid for all w(0) < 1 and all t > 0. 

Remark 5.1. The seemingly simpler modification (used by Le Jan and Sznit- 
man [8| for their uniqueness proof and by Bhattacharya et al. [3|) where a 
particle with label n produces two offspring each with label n—l, leads to 
the explicit iterative scheme u n = —u n +u„_ 1 . Unfortunately, the limit of 
u n (t) for large t, as n — > °°, fails to exist for w(0) < — 1. 

The semi-implicit approximation scheme works for other polynomial 
non-linearities. For example, if one considers u = — u — m 3 , the approxima- 
tion scheme u n = —u n — u^ ^Un, where each particle with label n branches 
into three particles, one with label n and two with label n— 1, is convergent 
to the true global solution for any initial condition. 

5.2. A general approximation scheme. The aim is to define a sequence 
of approximations Xk (0 t0 our abstract system of ODEs (I2.ll) . These ap- 
proximations will have a stochastic representation without any integrability 
problems. 

Rather than construct a particle system with labelled particles as de- 
scribed in the previous section, we put the modification into the evalua- 
tion operators. We claim there exists a sequence of evaluation operators 
Rn,t '■ & — ► C r satisfying the following implicit relations on Mo,;] < °°: 

I otherwise 
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and, for n > \,R nt (T) equals 
(5.1) 
fXk (O) 

7ke ('-'») 
CfR n , t - t9 (n (T)) 



t<b > t, 

tQ < t, death at 0, 
t% < t, flip at 0, 



C h B kf) , k k {R n ,t-t t) {ni('r)),Rn-l,t-t 9 (M' r ))) h <t branch at 0. 



The existence of R n j can be established exactly as in Lemma 1331 The intu- 
itive link with the labelled particle picture in the last section is that R n j(T ) 
corresponds to the evaluation operator applied to the tree started at a particle 
with label n at position k. 



(i) 



(i) 



(0) 

1 



(0) 

1 



(2) 



(2) 



(1) 



(2) 



(1) 



(1) 



(1) 



(0) 



(0) 

1 



Figure 3. The tree on the left can be pruned by starting, 
top right, with a particle labelled ( 1 ) , or, bottom right, by a 
particle labelled (2). Circles mark the death via pruning. If a 
label larger than (2) is given to the starting particle, the tree 
is unpruned. 

The implicit relation implies that if Afy^] < m then R njt (T) = Rt(T) 
whenever/? > m. Moreover when R nt (T) ^R t (T) theni?„,( , r) =0. Thus 
there exist increasing sets Q, n t C so that 

(5.2) RnA (T k)=Rt('rk)la n>t and {N M < <*>} c (J*V 

n 

We now define the stochastic representation using these modified evalu- 
ations by 



(5.3) 



xi\t) = E k [R n , t ]. 



The fact that this expectation is always well defined is part of the following 
result. 

Proposition 5.2. Suppose that %(0) e £°°(C r ) and y E L°°([0, r],^°°(C r )). 
Then the expectations in (15.31) are well defined and X k {t) are the unique 
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L°°([0, T],£°°(C)) mild solution to the following approximation scheme 
(5.4) 

%k = — ^kXji 5 

l,meZ rf 

with initial condition X^'(O) = %k(0) for all k G Z d and nGN. 

Proof. The local existence and uniqueness of solutions for the approxima- 
tion scheme, follows from the same methods as in the proof of Theorem 
13.11 plus an inductive argument in n > 0. The fact that solutions are glob- 
ally defined follows, again by induction, from the simple estimate 

tf(0l<llx(0)||~ + sup 11711-+ 

te[0,T] 



+ l k | f e-^( f - i ')(C / + C,||%("- 1 »|| oo )||x("»|| oo ^, 



which, using induction and Gronwall's lemma, easily gives boundedness of 
||jC (n) ||oo in each interval [0,T]. 

In order to prove that the stochastic representation (15.31) is well defined, 
we use a comparison argument, as in Section |4| The comparison equation 
for the approximation scheme is given by 

and the evaluation Ek|i?„ f | is finite as long as the %k are finite. But this 
follows by the same arguments as in first part of this proof. Again E k \R n>t \ < 
Xk < C for all k G Z rf and all t G [0, T] with constant C depending only on 
7\x(0),andy. 

Finally, the expectations Ek[/?„ ;/ ] do form the unique solution to the ap- 
proximation scheme by conditioning on the first branch of the tree as in 
Theorem 13^61 □ 



In the integrable case, that is where Ek|i?r| < °°, we have immediately 
from d5~2b that 

limE k [R n , t ] = E k [R t ]. 

In particular, when the expectations E are bounded over t G [0, T] and 
k G Z d this implies the solutions of the approximation scheme converge 
to those of the original system (I2.ll) . Our interest, however, is in the non- 
integrable case and we aim to show that convergence of the approximation 
scheme directly and deduce that the limit lim„^c« Ek[R n ,t] exists and defines 
a stochastic representation for all times t > 0. 
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5.3. Global convergence of the stochastic approximation. The aim of 

this section is to give a few details of one example where the approxima- 
tion scheme defined by the pruned representation converges, even when the 
direct stochastic representation fails to be integrable. In contrast to the pre- 
vious section, we use PDE methods. The convergence depends crucially on 
the equation and how the pruning is done, as not all approximation schemes 
will converge globally. 

For simplicity we work with the one dimensional Burgers equation with 
forcing (12.51) . In Subsection 12.21 we recast the equation into our abstract 
form by considering the weighted Fourier coefficients 

Xk(f) = WkKk(f) 

where, as in section l2~2l the weights are given by vv k = (1 V |k| a ) for some 
a > 1. If we assume the Fourier coefficients of the initial condition satisfy 

(5.5) sup{| Mk (0)|w k } <oo 

k 

and the forcing function / satisfies 

(5.6) sup sup {|/kWkk}<°°, foralir>0, 

k te[0,T] 

then proposition 15 .21 implies there is a unique global solution % k gi ven 
by (15.31) . to the approximation equations (I5.4I) . 

Theorem 5.3. Assume, in addition to (15.51) and (15.61) . that w(0) G H 1 and 
f G L^ c ([0,°°),L°°). Consider the pruned approximation of the previous 
section. Then the limit 

(5-7) Xk(0 = lim x£ } (0 = Km E k [/?„,,] 

exists for allt>0 and all k G Z and defines a global solution of the Fourier- 
transformed Burgers equation. 

Proof. Define 

Since % (n) is bounded we may reconstruct from these coefficients the func- 
tion 

« w (0=E«4" , (0e lfa . 
kez 

Using the representation of Proposition 15.21 we see that, on the level of 
PDEs, M (n) solves the approximation scheme given by 

d t u {n) = d 2 x u (n) + d x u {n) u { "- x) + /, 
«W(0) = k(0). 
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Fix T > and set C(f, T) = sup fe j r i ||/(?) We first use a maximum 
principle argument to show 

(5.8) sup \\u (n \t)\\ L ~ < ||«(0)|| L - + C(/,r)r forallneN. 

fe[0,r] 

We now derive an a priori estimate for the solution. The following calcu- 
lation applies to sufficiently smooth functions and standard approximation 
techniques imply that the resulting bound holds for the solutions above. 
Using (I5.8I) . we find 



fill 

d x u in) \\ 2 L2 = -\\%uM\\l 2 -J d x u ln) d n - l) d 2 x u (n) dx 

- / d 2 x u (n) fdx 
Jo 

< - \\d 2 x U {n) f L2 + C\\d 2 x U ( " ] ||^ 2 + C\\d 2 x U (n) \\ L 2 



where we have used the Poincare and Cauchy-Schwartz inequalities. Note 
that the constant C > depends only on T, C(f, T), and w(0). Thus we find 
another constant, also denoted C, such that for all n e N 

rT 

sup \\u in) (t) \\h <C, / \\u [n) {t)\\ 2 H 2dt<C 
ts[o,r\ ~ Jo 

and 

J \\d t u {n) {t)\\ 2 L 2dt<C. 

We now use standard methods to show that we have a solution of the limit- 
ing equation (cf. for example Temam [ 16 1). Indeed by compactness results, 
there is a subsequence («A;)feN> such that u nk — > u weakly in L 2 ([0, T],H 2 ) 
and H l ( [0, r],L 2 ), and strongly in LP ( [0, T] , L 2 ) for any p > 1 . Thus u is the 
weak solution of Burgers equation, i.e. it solves the PDE in L 2 ([0, T],L 2 ). 
As weak solutions of the Burgers equation are unique, we can neglect the 
subsequence, as any limiting point of u (,)) defines the same solution u. Fi- 
nally, the convergence is strong enough, in order to have all Fourier coef- 
ficients convergent. Thus for all k G Z the Fourier coefficients of u are 
given by 

Mk (0=lim u i ; ) (t) = \imw k l x^(t). 

□ 

Remark 5.4. We point out that the assumptions of the previous theorem are 
by no means optimal. We have used a simplified method of proof, in order 
to provide an example in a simple context. In particular the constraint on 
the initial condition can be relaxed. Furthermore, using regularisation prop- 
erties of the PDE, we can always get sufficiently smooth initial conditions, 
if we wait a small amount of time. 
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